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We argue that the smah fraction of neutrinos that undergo direction-changing scattering outside of 
the neutrinosphere could have significant influence on neutrino flavor transformation in core-collapse 
supernova environments. We show that the standard treatment for collective neutrino flavor trans- 
formation is adequate at late times, but could be inadequate in the crucial shock revival/explosion 
epoch of core-collapse supernovae, where the potentials that govern neutrino flavor evolution are 
affected by the scattered neutrinos. Taking account of this effect, and the way it couples to entropy 
and composition, will require a new paradigm in supernova modeling. 
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In this letter we point out a surprising feature of neu- 
trino flavor transformation in core-collapse supernovae. 
These supernovae have massive star progenitors which 
form cores which collapse to nuclear density and pro- 
duce proto-neutron stars. The gravitational binding en- 
ergy released, eventually some ~ 10 % of the rest mass 
of the neutron star, is emitted as neutrinos of all fla- 
vors in a time window of a few seconds. Diverting a 
small fraction of this neutrino energy into heating can 
drive revival of the stalled core bounce shock [lH7] creat- 
ing a supernova explosion and setting the conditions for 
the synthesis of heavy elements [H [31-15] • However, the 
way neutrinos interact in this environment depends on 
their flavors, necessitating calculations of neutrino flavor 
transformation. These calculations show that neutrino 
flavor transformation has a rich phenomenology, includ- 
ing collective oscillations [T0H38] . which can affect im- 
portant aspects of supernova physics [El [161 119^231 (271 - 
[^ [5T1 \T2[ l55Hi5] . For example, neutrino-heated heavy 
element r-process nucleosynthesis [11H15] and potentially 
supernova energy transport above the core and the ex- 
plosion itself [ni [23 [13 could be affected. 

All collective neutrino flavor transformation calcula- 
tions employ the "Neutrino Bulb" model, where neutrino 
emission is sourced from a "neutrinosphere" , taken to be 
a hard spherical shell from which neutrinos freely stream. 
This seems like a reasonable approximation because well 
above the neutrinosphere scattered neutrinos comprise 
only a relatively small fraction of the overall neutrino 
number density. However, this optically thin "halo" of 
scattered neutrinos nonetheless may influence the way 
flavor transformation proceeds. This result stems from a 
combination of the geometry of supernova neutrino emis- 
sion, as depicted in Fig. [l| and the neutrino intersection 
angle dependence of neutrino-neutrino coupling. 

Neutrinos are emitted in all directions from a neutri- 
nosphere of radius R^, but those that arrive at a loca- 
tion at radius r, and suffer only forward scattering, will 
be confined to a narrow cone of directions (dashed lines 
in Fig. [I]) when r 3> i?,y. In contrast, a neutrino which 
suffers one or more direction-changing scattering events 



FIG. 1: Supernova neutrino emission geometry. 

could arrive at the same location via a trajectory that 
lies well outside this cone. 

Following neutrino flavor evolution in the presence of 
scattering, in general, requires a solution of the quan- 
tum kinetic equations (50H52J . However, the rare na- 
ture of the scattering that generates the halo suggests 
a separation between the scattering-induced and coher- 
ent aspects of neutrino flavor evolution. In the coherent 
limit the neutrino-neutrino Hamiltonian, H^^, couples 
the flavor histories for neutrinos on intersecting trajec- 
tories [Ml HH [Sni [53]. As shown in Fig. [T] a neutrino 
Vi leaving the neutrinosphere will experience a potential 
given by a sum over neutrinos and antincutrinos located 
at the same point as neutrino vi: 

Hyu = V2Gf ^ (1 - C0s6'ia)n^,a \i^u,s) (V'l^.al 
a 

- V2Gf ^(l-C0s6lia)rip,a|V'P,a) (-iAp.al, (1) 
a 

where the flavor state of neutrino Vg, is represented by 
\ipv,a), and 0ia is the angle of intersection between Vi 
and neutrino or antineutrino v^/i^a.- Here n^^a is the lo- 
cal number density of neutrinos in state a, and the 1 — 
cos 9ia factor disfavors small intersection angles, thereby 
suppressing the potential contribution of the forward- 
scattered-only neutrinos [TOl [TT] . Direction-altered scat- 



tered neutrinos may have larger intersection angles as 
shown in Fig. [ll and therefore can contribute signif- 
icantly to the Savor-changing potentials, despite their 
small numbers. 

In the mean-field, coherent approximation, neutrino 
flavor evolution is governed by a Schrodinger-like equa- 
tion [54], id\ip^^i)/dt = H\tp^^i), where t is an Affine pa- 
rameter along neutrino v^s world line, and H — Hy + 
He + H^i, is the appropriate neutrino propagation Hamil- 
tonian, with vacuum and matter components Hy and H,^, 
respectively, i?^^ can be split into two pieces: H^^^^, 
contributed by neutrinos (index j in Fig. II]) which propa- 
gate directly (straight lines) from the surface of the neu- 
trinosphere; and i?^^'°, contributed by neutrinos that 
suffer direction-changing scattering outside the neutri- 
nosphere (index k in Fig. fl]) and propagate coherently 
thereafter. To wit, H^^ = TJ^uib _^ ^haio^ 

The operators H})^^° and H^^^^ depend on the com- 
plex phases of the neutrino flavor states which contribute 
to them, so that the relative leverage of these opera- 
tors in determining flavor transformation at any point re- 
quires numerical calculations. Some conditions have been 
shown to give phase locking, while other conditions give 
phase decoherence [TH [Ml |331 I3H1 SI] • For the purpose of 
evaluating the validity of the Neutrino Bulb model, we ig- 
nore path length difference-induced phase averaging |44j 
and compute the maximum magnitude of the diagonal 
Hamiltonian elements, which we denote with |if^^'°| and 
|/f^""^|. A necessary condition for the validity of the 
Neutrino Bulb model is that \H^^^^\ > \H^!;^°\. 

A simple argument can be made about which vari- 
eties of spherically symmetric density profiles could ren- 
der the Neutrino Bulb model inadequate. Consider a 
series of spherical shells of matter stacked around the 
neutrinosphere. These shells are taken to isotropically 
scatter neutrinos, and, as we discuss below, neutral cur- 
rent neutrino-nucleon/nucleus scattering does just this. 
Some of these neutrinos will contribute number density 
and flavor information to the sum in Eq. [l] adding to 
H^^^^°. For a point r well outside of these shells, the 
number density of neutrinos being scattered to this loca- 
tion from a shell at radius r', multiplied by the average 
value of (1 — cos 0ik) for neutrinos coming from this shell, 
is ex p{r')6r' (r'/Ri,) , where 5r' is the thickness of the 
shell. As r' approaches r, the contribution from these 
shells is regulated. For small r', the shell contributions 
are regulated by the neutrinosphere. 

These considerations imply that when p(r') ex r' , 
the potential contributed by a given shell will be ex 
6r' /r' (X (Slog(r'). Any selection of logarithmically 
spaced shells with r' < r will cause each shell to make 
an equal contribution of neutrino number density at r. 
Physically, we might expect density features of size r' at 
radius r'. This matter density configuration will cause 
the ratio of | H^^^° \ / \ H^^^^ \ to remain fixed with increas- 
ing radius. 

To compare the contributions of the halo shells to 
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FIG. 2: Solid lines show matter density profiles and dashed 
lines the corresponding Neutrino Bulb (1%) safety criteria 
from Eq. [2] Black lines are for the late-time neutrino driven 
wind environment [15], green lines the neutronization burst 
0-Ne-Mg core-collapse environment [551156] . and red lines the 
Fe-core-collapse shock revival environment |57) . 



l-f^i'?""^|j we observe that the neutrinosphere (more pre- 
cisely, the transport sphere [35]) can be treated in the 
same spirit. The transport sphere is characterized by 
the neutrino optical depth, r, equal to unity. Requiring 
that the logarithmic shells above the neutrinosphere con- 
tribute much less than the neutrinosphere itself results in 



p (r) < pr=i 



Ru 



(2) 



Early in the explosion epoch, the transport sphere corre- 
sponds to physical radii R^ ~ 30 — 60 km and densities 
Pr=i ~ 10'" *° ^^^ gcm-3 [58]. In fact, Eq. [2] assumes 
that the thickness Ai?,^ of the neutrinosphere is ^ i?^, 
whereas models show that /S.R1, < R^^ implying a more 
stringent constraint by a factor of AR^/R^. 

If |iJ^^'°|/|i/^;^"^| < 1% is taken as the hmit where 
ijj,'^'° can be neglected, then by Eq. ^ the range of den- 
sity profiles for which the Neutrino Bulb model is likely 
to be adequate is p (r) < 0.01 x Pt=i {R„/r) . As long 
as the matter density in the supernova remains below 
this limit, there is no danger that the fractional poten- 
tial contribution from the scattered halo, |iJj'^^'°|, will 
grow above 1 %. Fig.[2]shows the density profiles for sev- 
eral core-collapse supernova environments alongside the 
corresponding 1 % safety criterion for each profile. 

As is evident in Fig. [2| the Fe-core-collapse shock re- 
vival environment will have a significant scattered halo. 
Even though the 0-Ne-Mg core-collapse density pro- 
file [S5, 56J drops into the safe zone at r > 1000 km, these 
models nevertheless will have a significant scattered halo 
originating from shells at lower radius where the density 
curve is above the 1 % safety margin. Only late-time 
neutrino-driven wind models avoid scattered halo com- 
plication [TCll^ l^g i lM l lSgHST]. Fe-core-collapse models 
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FIG. 3: Left: Color scale indicates the density within the shock front in a ISM© progenitor core-collapse supernova 500 ms 
after core bounce, during the shock revival epoch [53- Right: Effect of the scattered neutrino halo for the matter distribution 
at Left. Color scale indicates the ratio of the sum of the maximum (no phase averaging) magnitudes of the constituents of the 



neutrino- neutrino Hamiltonian, \H,, 



I 1 1/^ 



to the contribution from the neutrinosphere \H, 



bulb I 



(e.g., the red curve in Fig. ^, in general, exhibit an av- 
erage density profile that is ex 



-(2 to 3) 



which means 



that l-ffi^^'"!/!^]^""^! is expected to increase with radius. 
Note, however, that though the relative contribution of 
the halo may grow with radius, at sufficiently large dis- 
tance from the proto-neutron star the neutrino-neutrino 
potential ceases to be physically important. 

Matter inhomogeneity, an essential feature of super- 
nova explosion models HHZl HTJ HH [33] , adds complexity 
to this issue. To study this effect we use the 2D mat- 
ter density distribution. Fig. [Sj taken from a supernova 
model derived from a 15 AIq progenitor [57]. This snap- 
shot corresponds to 500 ms after core bounce, during the 
shock revival epoch, after the onset of the SASI [H [S] . 
We mock up a full 3D density profile by cloning the 2D 
profile into a 3D data cube. Starting with an initial flux 
of neutrinos from the neutrinosphere [64] . and taking all 
baryons to be free nucleons, we use the full energy de- 
pendent neutral current neutrino-nucleon scattering cross 
sections |65| to calculate the number flux of neutrinos 
scattered out of each spatial zone and into every other 
spatial zone (retaining the necessary information about 
relative neutrino trajectories between zones). We com- 
pute the magnitude of |-ffj^^'°| at each location in the 2D 
slice that comprises the original density distribution. 

In this example calculation the scattered halo is taken 
to be composed of neutrinos which have suffered only a 
single direction-changing scattering. Because the halo re- 



gion is optically thin for neutrinos, multiple scatterings 
become increasingly rare with radius and do not have a 
geometric advantage in their contribution to |-ff^^'°| rel- 
ative to singly-scattered neutrinos. Neutrinos which ex- 
perience direction-changing scattering that takes them 
into the same cone of directions as neutrinos forward 
scattering from the neutrinosphere are counted as con- 
tributing to the halo (these neutrinos contribute ~ 10~^ 
of the halo potential). As before, we neglect the effects 
of neutrino flavor oscillations. Fig. ]3] shows the results 
of this calculation out to a radius of r = 2000 km. Dis- 
turbingly, neutrinos from the scattered halo in this 2D 
model nowhere contribute a maximum magnitude less 
than 14 % of the neutrino-neutrino potential magnitude, 
and in many places contribute 90 % or more of the total. 
Fig. ]3] shows that matter inhomogeneities generate large 
corresponding scattered halo inhomogeneities. 

The inhomogeneity of the scattered halo is increased 
by several scattering processes which have been omitted 
from this illustrative calculation. We did not include 
neutrino-electron scattering. This scattering process has 
smaller cross sections and relatively forward peaked an- 
gular distributions and therefore produces a subdominant 
contribution to |-ff^j^'°|. What is more important is that 
our calculation leaves out what is likely the dominant 
source of neutrino direction-changing scattering in the 
low entropy regions of the supernova envelope: coherent 
neutrino-nucleus neutral current scattering. 



The cross sections for this process scale as the square 
of the neutrino energy and square of the nuclear mass 
number A. In fact, since the proper number density of 
nuclear targets is oc A~^, but the coherent scattering 
cross section c>c A^, the overall scattered halo potential 
contribution stemming from this process is oc A. This 
process, like neutral current neutrino-nucleon scattering, 
is flavor independent and flavor preserving, simply chang- 
ing neutrino direction. 

Since the heavy nucleus mass fraction, and the dis- 
tribution of nuclear mass numbers, can depend sensi- 
tively on the entropy and electron fraction [66] . coherent 
neutral current scattering could couple neutrino flavor 
transformation to macroscopic, multi-dimensional struc- 
tures in the supernova envelope. For example, the model 
shown in Fig. [3] has relatively lower entropy, downward- 
flowing, higher nuclear mass fraction matter; and higher 
entropy, upward-flowing plumes, with relatively lower nu- 
clear mass fraction. This could produce a scattered halo 
with a complicated 3D geometry and flavor content, cre- 
ating a non-trivial enhancement to the inhomogeneities 
evident in the scattered halo potential shown in Fig. [3j 

Because scattering processes are energy dependent, 
neutrinos in the scattered halo possess different energy 
spectra than forward-scattered only neutrinos. Fur- 
thermore, the flavor content of the scattered halo will 
not match that of neutrinos emerging from the ncutri- 
nosphere. For example, consider emergent Vc and Dc 
number fluxes that are equal, yet have different energy 
spectra. Taken alone, unmolested by neutrino flavor os- 
cillations, these fluxes give l^^'j'""^] = 0. However, passing 
through the energy-dependent scattering processes, they 
yield |_ffh^'°j ^ 0. 

Anticipating the course of neutrino flavor evolution 
in this environment is clearly challenging. \Hc\ is by 
far the largest contribution to \H\ during the shock 
revival epoch. However, considering the Bulb neutri- 
nos alone, current coherent calculations show that neu- 
trino collective flavor oscillations can proceed despite a 
large matter potential [13 [HI [13 HO]- The criterion 
for the matter suppression of collective oscillations [67] , 
A I He I ^ A|i?v|j where A denotes the dispersion in these 
potentials for Bulb neutrinos, is not met where the mat- 
ter densities in Fig. ^drop below ~ 10^ *° ^ gcm~'^ (the 
wide range is due to the geometric dependence of A|i?c|)- 
In a further complication, neutrinos from the spatially 
extended scattered halo could arrive at a given loca- 



tion along many different trajectories with different path 
lengths, so that significant neutrino oscillation phase av- 
eraging [H] could come into play. This has been shown 
to suppress collective oscillations in some conditions [33] . 
Inhomogeneity and the intersection angle dependence of 
the neutrino-neutrino interaction may make phase aver- 
aging incomplete. Ascertaining the role of decoherence 
and phase averaging processes requires detailed calcula- 
tion with specific supernova models [SSj. Even if col- 
lective oscillations are found to be suppressed at small 
radius, they may be operating, e.g. above the shock, 
because the halo extends the collective oscillation region. 

Though validating coherent flavor transformation 
studies for late-time neutrino driven wind models, our 
calculations demonstrate the potential inadequacy of 
these treatments in an environment important for the un- 
derstanding of the supernova explosion mechanism and 
nucleosynthesis. Ultimately, the scattered halo changes 
the nature of the neutrino flavor transformation prob- 
lem: it broadens the region influencing flavor evolution 
from just the neutrinosphere to a much larger fraction 
of the supernova envelope; and it introduces essential 
multi-dimensional effects. The standard Neutrino Bulb 
model by its nature is an initial value problem at each 
radius r, while the scattered halo makes it necessary to 
consider how flavor transformations at large radii can 
feed back into the evolution at smaller radii. With this 
additional source of nonlinearity, qualitatively new phe- 
nomena could, in principle, occur. Further, the extended 
scattered halo can couple neutrino flavor evolution to 
the nuclear composition and complex 3D flow geome- 
tries which are characteristics of the supernova explo- 
sion epoch. A self-consistent solution of this problem 
likely will demand new computational capabilities and 
approaches. Given the importance of neutrinos and the 
supernova phenomenon for so many aspects of our under- 
standing of the cosmos, it may be that there is no choice 
but to seek such a solution. 
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